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We discuss the mean-field theories obtained from the leading order in a large-TV approximation for one- and 
two- component dilute Bose gases. For a one-component Bose gas this approximation has the following prop- 
erties: the Bose-Einstein condensation (BEC) phase transition is second order but the critical temperature Tc 
is not shifted from the non-interacting gas value To. The spectrum of excitations in the BEC phase resembles 
the Bogoliubov dispersion with the usual coupling constant replaced by the running coupling constant which 
depends on both temperature and momentum. We then study two-component Bose gases with both inter- and 
intra- species interactions and focus on the stability of the mixture state above Tc- Our mean-field approxima- 
tion predicts an instability from the mixture state to a phase-separated state when the ratio of the inter-species 
interaction strength to the intra-species interaction strength (assuming equal strength for both species) exceeds a 
critical value. At high temperature this is a structural transition and the global translational symmetry is broken. 
Our work complements previous studies on the instability of the mixture phase in the presence of BEC. 

PACS numbers: 67.85.Fg,03.75.Hh,05.30.Jp 



I. INTRODUCTION 

With the capability of tuning inter-particle interactions and 
the realization of superfluid phases (see U] |2| for reviews), 
cold atoms offer a window of unparalleled promises onto 
many-body physics. While the cold atom prospect of study- 
ing quantum criticality has attracted much attention (T], the 
finite temperature thermodynamics |4| offers an equally fer- 
tile ground for explorations of fundamental questions. The 
power of effective field theory has also proven to be useful in 
treating other aspects of many-body physics in cold atoms [Sj. 
Recently developed technologies for accurate temperature de- 
termination |6| and for creating stable, flat (bulk-like) trap- 
ping potentials |7 | and ring-shape potentials |8 | provide some 
of the necessary tools for probing thermodynamics and phase 
transitions of cold atoms. Phase separation, the demixing phe- 
nomenon that spontaneously breaks the global translational 
symmetry, was the second transition after Bose-Einstein con- 
densation (BEC) to be observed in the cold-atom laboratory 
im. The phase separation transition that provides a paradigm 
of second-order scaling physics in ordinary finite temperature 
phase transitions was observed in a mixture of dilute BECs 
near zero temperature. With cold atom technology, the dy- 
namics of a zero temperature miscible-immiscible transition 
of bosonic superfluid mixtures, first discussed in studies mo- 
tivated by the plans of creating liquid ''He-^He mixtures | lOj, 
can now be studied in trapped atoms ifTTll . The magnetically 
controlled Feshbach resonance ||2] [l2l \T3l provides a direct 
trigger and promises a useful probe of the interaction depen- 
dence of the phases and phase boundaries. Here we describe 
the finite temperature phase separation transition above the 
critical temperature Tc of BEC in a two-component boson 
mixture. 
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The thermodynamic descriptions of cold atoms encoun- 
tered fundamental challenges posed by the inherent limita- 
tions of the gapless and conserving approximations lfT4l of 
interacting bosons. As a consequence, the treatments of the 
BEC-transition in a single component boson gas generally ob- 
tain a transition that is first order whereas it is known to be 
second order [4]. The long-standing problem of the interac- 
tion dependence of the Tc of a single component BEC (see 
|4| for a review) was then discussed in the low density limit 
by intricate reasoning tailored to the computation of Tc (only) 
ifTSl [T6l . More recently ifTTl [TSl . we developed an auxiliary 
field description that can overcome the obstacles of the con- 
serving/gapless approximations and reproduce the correct or- 
der of the BEC-transition at the mean-field level. This for- 
malism, the Leading Order Auxiliary Field (LOAF) approxi- 
mation, introduces two composite fields, one to describe the 
density 0*0 and one to describe the anomalous density (fxp, 
where denotes the bosonic field. This treatment also pre- 
dicted a low density limit of the T^-dependence upon the scat- 
tering length consistent with previous work based on large-A^ 
expansions |17| while providing a complete description from 
which all thermodynamic quantities can be computed 1 1 8 1 . 
Above the critical temperature, where the anomalous density 
vanishes, this description becomes identical to the usual large- 
N approximation ||T9| . Below, we describe a two-component 
Bose gas in the large-N approximation and investigate the sta- 
bility of the homogeneous (mixture) phase above the BEC 
transition temperature to the leading order To provide context 
and to gauge the performance of this method, we also derive 
and discuss the large-N predictions for the thermodynamics of 
a single-component finite-temperature Bose gas. 

Although the relativistic 0{N) model for self-interacting 
bosons has been the subject of many papers ||2QH23 I. the non- 
relativistic version which introduces N replicas of the 0{2) 
or equivalently U{1) symmetry of the non-relativistic dilute 
gas effective theory has not been discussed in great detail in 
the literature. The main use of the large-N expansion in BEC 
theories has been to discuss a large-N approximation near the 
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critical temperature Tc as done by Baym et. al. fTSl and 
Arnold and Tomasik li24J as well as in the work of Braaten and 
Radescu ifTSl and has been reviewed by Ref. [4]. In this pa- 
per we will discuss the broken symmetry aspects of this prob- 
lem following the classic paper of Ref |20|. For pedagogical 
purposes we will use in this paper a more general method of 
introducing auxiliary fields discussed by Refs. |[T9l |25l and 
briefly discuss the (equivalent) Hubbard-Stratonovich method 
that we used in ifTTHTSl and which was also used in 1201 . The 
Hubbard-Stratonovich method relies on one having quartic in- 
teractions so it may be appropriate to introduce the atomic 
physics community to the more general method of introducing 
auxiliary fields that can be used for arbitrary polynomial (as 
well as non-polynomial) interactions that preserve the replica- 
tion symmetry |26|. In this paper we will first derive the large- 
N expansion for a single component Rose gas. The large- 
expansion when evaluated (as we do here) at iV = 1 is equiv- 
alent to choosing 9 = Q in the auxiliary-field approach dis- 
cussed in detail in Ref 1 18 1 so it does not include the anoma- 
lous density explicitly. In Refs. ifTTlfTSjl we introduced a loop 
counting parameter e, which is identical to the loop counting 
parameter above T^. where the anomalous density van- 
ishes. 

What we will show is that for a one-species gas of bosons, 
the leading order in our large-iV approximation leads to a non- 
perturbative (in coupling constant) mean-field theory with rea- 
sonable features. As in the more sophisticated LOAF treat- 
ment, the leading-order large-iV approximation also predicts 
the correct second order EEC transition. Importantly, we will 
show that the large-A^ theory does lead to a Bogoliubov-like 
spectrum for temperatures below Tc because of mixing be- 
tween the fluctuations of the boson and the composite field 
when there is a broken symmetry. However, a shortcoming of 
this expansion is that it does not predict in the leading order 
in a shift in the critical temperature from that of the free 
gas, a feature that is shared by the Popov approximation ifTSl . 
This defect is rectified at the mean-field level by also includ- 
ing an auxiliary field for the anomalous condensate as in the 
LOAF approximation. Since the LOAF approximation leads 
to the same result as the large-A^ approximation above Tc and 
we are interested in the stability of the mixture phase of a 
two-component Bose gas above Tc, we will study the simpler 
large-A^ approximation here, which ignores the contributions 
from the anomalous density in the condensate regime in the 
leading order As summarized in Ref |4| once the cor- 
rections (and higher order corrections) to the self energy of 
the boson propagator are included, one does find a shift in Tc 
so that the expansion at higher orders does include the 
effects of the anomalous density. Finally, we note that if com- 
parisons with experimental results in an inhomogeneous trap 
are needed, one may use the local density approximation 1271 
to include possible inhomogeneity effects. 

This paper is organized as the following. Section|ll]presents 
our large- approximation for a single-component interacting 
Bose gas both below and above the BEC transition tempera- 
ture. The excitation spectrum in the BEC phase will be ana- 



lyzed in details. Section III shows our large-iV approximation 
for two-component Bose gases in the mixture state as well as 



the phase-separated state. A phase transition between the two 
states is found in the normal phase and we present a phase 



diagram from our theory. Section IV concludes our work. 



II. 



LARGE-iV THEORY FOR A SINGLE-COMPONENT 
BOSE GAS 



The partition function of a single-component Bose gas can 
be given a many-body theory path-integral representation 14] 



(1) 



where we are using the Matsubara imaginary time formalism. 
/3 = /i is the chemical potential, and V is the vol- 

ume of the system. The Euclidian action S[(j), 4>*; V, /i, /3] is 
given by 



where we have introduced the notation 



dr. 



(2) 



(3) 



For a dilute Bose gas the effective field theory for the problem 
can be describe by the Euclidian Lagrangian density IH 



£ = - [0*9,^- 



2m 



2m 



(4) 



Here A is the 
will discuss its 
grangian density 

H = jd^x 



bare coupling constant and we 
renormalization later This La- 
corresponds to the Hamiltonian 



+ 



In 



f \X{4>*4>) 

To determine the 



27n ^ 2m ' 

what follows we set ^ = 1 and = I 
finite-temperature effective potential for the theory we will 
be interested in the generating functional for the connected 
correlation functions, lnZ[j], where 



(5) 

The large-A^ expansion is a combinatoric trick that re- 
organizes the Feynman diagrams of the theory in a non- 
perturbative fashion. First it sums the loops (bubbles) con- 
tributing to the scattering amplitude. Calhng this bubble sum 
a "composite-field propagator", one then re-sums the theory 
by implementing a loop expansion in terms of the number of 
composite-field propagator loops in a diagram. This way of 
organizing the Feynman diagrams of the theory can be ac- 
complished formally by introducing N copies of the original 
theory which is equivalent to introducing a "color" index with 
N components, or in some cases extending a theory with a 



3 



0(1) or 0{2) symmetry to an 0{N) symmetry. For the La- 
grangian density (H]) one then introduces (see below) a com- 
posite field oe = ^ 4'n't'n into the theory by introducing 
formally a "1" into the generating functional Z[j] shown in 
Eq. (|5]l by using a functional expression for the delta function 
enforcing the definition of a. Equivalently the composite field 
can be introduced using a Hubbard-Stratonovitch transforma- 
tion as shown in |17 18 , 20|. As shown below this converts 
the quartic self interaction into a trilinear interaction that is 
quadratic in the original (j) field. This allows one to perform 
the path integrals over the original fields 0„ in the generating 
functional Z[j] exactly while keeping a fixed. At this stage 
one could then determine a by using its definition and obtain 
the Weiss self-consistent mean-field theory ||29]| . The beauty 
of the path integral approach is that this mean-field theory is 
the first term in a complete resummation of the theory in terms 
of loops of higher and higher numbers of the composite-field 
propagators flT]. Having N copies of the original theory or 
extending 0{2) to 0{N) introduces a small parameter 
into the theory which allows one to perform the remaining 
path integration over the composite field a, which arises from 
inserting the formal expression for the delta function into the 
path integral, using Laplace's method (or the method of steep- 
est descent). 

To make this procedure explicit, we first make N copies of 
the original field Il20ll28l by generalizing (j) 0i, • • • ,(t>N, 
rescale the coupling constant A — > X/N, and define <i> — 
{(j>i ,(f>l,--- , (pN , 4'*n)'^ ^'^'^ external sources J so that the 
generating functional for the correlation functions is given by 



N 



(6) 



where the action S is given by 



\dx\ 




(7) 

Here J = ( ji , ]{ , j2 , J2 , ' ■ ' , ^iliV is the source coupled 
to Go ^ = diag{h^+\h^~\ ■■■ , {N identical 

copies), Gq ^ = G(7^ — diag{ii, /i, • • • , /i, /i) is the bare (non- 
interacting) Green's function, and /i'^^ = ±9^ — V^/(2m). 
The classical value of the n-th field is ipn^c = {^/Z){5Z/6j'^). 
Details of the large-N approach and its applications to other 
fields can be found in Refs. 1 19 30, 31 1 

We then introduce the auxiliary field a = ^ 
to facilitate our resummation scheme outlined above by in- 
serting the following identity inside the path integral for the 
generating functional ([TJ using a formal integral representa- 
tion of the Dirac delta function 



1 = J 'Da5{a — 



A 



N 

N ^ 

n=l 



N 



X{a- 



N 



N 



(8) 



Here = l/(27ri) is a normalization factor and the x in- 
tegration contour runs parallel to the imaginary axis as dis- 
cussed in Ref. |fT9l . This representation allows one to replace 



"Y^^i^i ^n4>n by {N/\)a in S inside the path integral. Let 
Go ^ = + diag{x, X, ' • • , X, x)- It is now possible to 
perform the quadratic integral over (pn exactly to obtain a new 
effective action that (because of the large factor A^) can be 
evaluated by Laplace's method. After integrating out 0„, one 
has 



Z[J,S,K] 



VxDae 



^ff 



(9) 



where we have added sources for the auxiliary fields x and a 
and 



[dx] 



1 . TV 



N , N 



-TrlnG^^ - (Sx + Ka) 



(10) 



The Tr In G^^ term comes from the Gaussian integration over 
the bosonic fields (see Eq. (fTTji). Note that the first term and 
the Tr In Gq'^ term are just A^ copies of the U{1) theory so 
they are of order A^. We can also rescale the sources S and 
K to be proportional to A^ so that a large parameter A^ is in 
front of the entire action. This enables us to evaluate the re- 
maining integrals over x and a by Laplace's method (or the 
stationary phase approximation). The resulting expansion is a 
loop expansion in the composite field propagators for x and a 

ca. 

The leading order in large-A^ expansion is obtained by just 
keeping the contribution to Z evaluated at the minimum of the 
effective action Seff (i.e. the stationary phase contribution) 
lEOlllII. Note that 5S'e///(5j* = -0„,c- Although we are 
interested in the theory with A^ = 1, for many problems the 
large-N expansion (which is an asymptotic expansion) gives 
qualitatively good results at A^ = 1 at leading order and the 
corrections at next order bring one closer to the exact answer 
This was seen in the calculation of Tc (in a slightly different 
context) in Ref. Il24l . 

Exactly the same large-N expansion can be obtained from 
completing the square in a shifted Gaussian integral. This is 
the well known Hubbard-Stratonovich transformation which 
is useful when the interactions are only quartic in nature and 
is based on the identity (here given for multi-dimensional in- 
tegrals) 



dxidx2 ■ ■ ■ dXr. 
(27r)»/2 

[detM]-i/2 



■ exp 



exp 



(11) 



I integral becomes 



On a lattice, with the substitutions ji - 
Sij/X, Xi — > a{i) we find that the 
quadratic, but we now have to be able to perform the result- 
ing (path) integration over the composite field a{x), which is 
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again done by the stationary-phase approximation. The result- 
ing large-N expanded effective action is the same as one ob- 
tains using the more general method of introducing the com- 
posite field a once one eliminates the Lagrange multiplier 
field x(x) from the problem, as will be shown below. 

From the Legendre transform of S^f / one obtains the gen- 
erating functional of the one particle irreducible graphs, which 
is the grand potential r[0, x, a] 1 2 8j i31j ,32J . Explicitly, 



r = / [da;](J^$c + Sxc + Ka,) + S, 



(12) 



where now <i>c,Xc,ac stand for the expectation values of 
<&, X, «■ We define the effective potential for static homo- 
geneous fields 0„,x,a as Ve// = YjNV^. Note that the 
Legendre transformation introduces the expectation values of 
0„ and 0* in F and T4/ / via SV ^ — in or, equivalently, 
J ~ Gq ^$ for the expectation values. Now that we have ob- 
tained the leading order approximation, we will set iV = 1 so 
that we are addressing the real dilute gas which has an C/(l) 
symmetry. At the leading order we find : 



2^ A^"+l^"'-X^" + 

1 



(13) 



Here we dropped the subscript c for the expectation value, 
$ = (0,0*)^ for the iV = 1 case, and G(7"^[x] has been 
reduced to a 2 x 2 matrix which depends on x- The next order 
in the 1 jN expansion involves the Gaussian fluctuations in the 
auxiliary fields a, x and will not be included here. 

The broken-symmetry condition is determined from the 
condition that we have found the true minimum of the ef- 
fective potential: <5Ve/ / — 0, which becomes x0 = 0. 
This imposes the following conditions: (1) In the normal 
phase = and a finite x is allowed and (2) in the broken- 
symmetry phase, (\) is finite so x = 0. 

After Fourier transforming and summing over the Matsub- 
ara frequencies, the last term of 14// becomes X^fcl'^fc/^ + 
(l/^)ln(l-e~'^"'')], where Wfc = efc + xande^ = k'^/{2ra). 
One can eliminate the Lagrange multiplier field x by using the 
minimum condition 5Vej f /5a ^i). Explicitly, 



X 



-/i + a. 



(14) 



Then one obtains 

K// = 



k 



2A 



a 



)]• 



(15) 



Here = efc — /i + a. 

We need to renormalize the theory because Eq. ([TSj is ultra- 
violet divergent. The renormalized coupling constants can be 
defined from the effective potential and it is the value of the 
scattering amplitude at zero energy- and momentum- trans- 
fer or equivalently l/Xu = S'^V^f j /SaSa = 1/A+ (finite 
polarization terms). The polarization terms can be shown to 



vanish at zero temperature, so if we define the renormalized 
coupling constant at T = then A/j — A, and one may write 
A = inh'^a/m, where a is the s-wave scattering length at zero 
temperature. 

The renormalization of the chemical potential fi can be seen 
more clearly if we rewrite Vi.efS in terms of x for the mo- 
ment. The unrenormalized effective potential is given by 



eff 



Vo + X't>*'f> 

E 



efc + X 
2 



1 

2A 
1 



ln(l 



)]■ 



(16) 



Here I^) is the unrenormalized vacuum energy. In the clas- 
sical theory —dVef f / dx ~ m/-^- So defining hr/X for the 
quantum theory via 



eff 



dx 



A 



(17) 



and only keeping the infinite contributions from the quantum 
fluctuations, we find 



A 



A 



(18) 



Finally there are infinite contributions to the potential that 
are independent of the field values. These do not contribute to 
the equations of motion but can be rendered finite by defining 
a finite constant Vr q via 



Vi 



R,0 



/4 

2A 



Vo 



2A 



k 



(19) 



With our choice of renormalized parameters we obtain for the 
renormalized effective potential 



VR.eff = Viifi + x4'*4' - 

sr^ 1 



ln(l 



1 

2A 



(20) 



It is often convenient to change variables and write everything 
in terms of ub- Similar to Eq. ( [T4] l we introduce a/; = x + 
fiR. The renormalized effective potential density becomes 



Veff = (-M + a) 



a 
2A 



(21) 



Here we drop the subscript R and the vacuum energy. = 
ek - fi + a. 

In the normal phase, = 0. The equations of state are 
derived from SVeff/Sa = and —SVeff/6^ = p. Exphcitly, 



a 
A 



k 



n(wfe), p = 2j'^(wfe). (22) 

fe 

Here ujk = £k — P + a and n{x) = [exp{(3x) — is 
the Bose distribution function. We define ko = p^'^ and 
use k^^ as the unit of length. The BEC transition temper- 
ature of a non-interacting Bose gas with density p is Tq = 
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2TTh^ p^/'^ /[^/^ {■i/2)kBm] and we use ksTo as our unit of 
energy. 

In the broken-symmetry phase, is finite and the condition 
SVeff/Scj)* = requires that fi ^ a. Therefore ujk = Cfc, 
which is the same as the dispersion of a non-interacting Bose 
gas. This imphes that the single auxihary field large- N the- 
ory used here does not lead to a shift in the critical tempera- 
ture Tc from that of a noninteracting gas. One may estimate 
the shift of the critical temperature by including higher-order 
terms (see Refs. |S] [12 and references therein). A more so- 
phisticated mean-field theory in the BEC phase, the LOAF 
theory which contains two auxiliary fields, has been studied 
in Refs. fTT T8 1 and does lead to a shift in Tc in the lead- 
ing order. Since the LOAF theory leads to the same result as 
the large-A^ approximation above Tc and we are interested in 
studying the stability of the mixture state of a two-component 
Bose gas above Tc, we will confine ourselves to the simpler 
leading order in the large- approximation which ignores the 
contributions from the anomalous density. However, as we 
shall see below, the large- theory does lead to a Bogoliubov- 
like spectrum below Tc- Thus it is a qualitatively reasonable 
approximation even below Tc- We next define pc = (f>*(l) 
as the condensate density in the broken-symmetry phase and 
consider a Bose gas of density p. Then 6Veff/Sa — and 
-6Veff/Sp = pgive 



a 
A 



(23) 



Since ujk = ^k, the second equation implies that Pc/p as a 
function of T is insensitive to A in this theory. 

Figure[T]shows pc, p, a, and the minimum of Veff as func- 
tions of T. One important feature is that the BEC transition 
temperature is fixed at Tq because the dispersion Uk is iden- 
tical to the dispersion of non-interaction bosons regardless of 
the interaction strength. The transition is second order be- 
cause Pc is continuous at Tc = Tq. The Vg// that is plotted 
is the value of the effective potential at the minimum where 
SVeff/Scj) = SVcff/Sa = 0. A negative value of Veff corre- 
sponds to positive pressure so the system should be mechan- 
ically stable. The BEC condition requires a = p, which can 
be verified below Tc- 

Once we have found the correct ground state of this theory, 
it is important to calculate the propagators in this ground state. 
This has been addressed in detail in the relativistic 0{N)- 
model in Ref. EOl and here we will follow a similar procedure 
in the broken-symmetry phase. In the broken symmetry phase 
(BEC phase) the and x propagators mix. To calculate the 
propagators in the broken symmetry phase one needs to invert 
the matrix inverse Green's function that is obtained from the 
effective action which is the generator of all one-particle ir- 
reducible graphs and which is the Legendre transform of the 
generating functional (— hi Z). The broken symmetry ground 
state is described by x = and (0) — J'p'c > 0, where pc is 



found from solving Eqs. (23 1 

The effective action whose static part leads to Eq. 
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Figure 1 . Single auxiliary field large- A' theory of a single-component 
Bose gas. (a) pd p vs. T. (b) p and a vs. T. (c) Veff (evaluated at 
the minimum) as a function of T. Red and black curves correspond 
tofcoa = 0.04 and 0.1. 



given by 



[dx] 



1$ 



N{x + Pf 
2A 



^Tr InGn 



(24) 



Here Gi,^[x] 



iuJn + i^k.i^n + Wfc, • • • , -iiOn + 

= Cfc + X- Since we are interested 



= diag( 
LOk,ii^n + ^k) and ujk 
in the A^ = 1 case we will now set A^ = 1 and confine 
ourselves to the actual one-component Bose gas which has 
a [/(I), or equivalently 0(2), symmetry. Let us first look 
at the U{1) approach. Here we can use the U{1) symme- 
try of the theory to choose the vacuum expectation value of 
(/) to be real. Thus in the broken symmetry phase we let 
— \fPc + 0, 0* — ^fp~c + 0* . The term x4'* 4> becomes 
x(Pc + \JlVc<^ + \fPc(\>* + 00*). The inverse propagator in 
the (x, 0) sector is not diagonal because of the condensate pc- 
Let 5" = (0, 0*, x)^ ■ Then the fluctuations can be written as 
vj/tp-iyj/ The inverse propagator matrix 'D~^ is obtained by 
taking the second derivatives of the effective action and then 
evaluating these in the broken symmetry ground state where 
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X = and {(f)) 



5r, 



\ 2 ,50*(5x 2 (50(5x 



5d,54> 

1 ^Tel/ 

2 5^'(5i5 

1 ■sr,^.,^ 



1 ST^ff \ 

2 ' 



(25) 



The upper 2x2 submatrix is just (l/2)(iza5(— za;„+Wfc, iw,, 



'5x(a;)'5x(y) 



6(x — y) 
A 



1 



TrGo{x,y)Go{y,x) 



(26) 



Here x, y denote the imaginary time and spatial coordinates, 
6{x — y) is the four-dimensional Dirac delta function, and we 
have used SGq/Sx — ~ J Go{SGq^ /Sx)Ga- After a Fourier 
transform this becomes 



SxSx 



(K) 



---B{K,T). 



(27) 



HereB(i^,r) = (1/2) TrGo(Q)Go(A' + Q), if - 

(iw„,fc), Q = iin^,q), J2q = (l//5)E,.Eq with w„ and 
ilj, being bosonic Matsubara frequencies. The expression for 
B{K, T) will be given shortly. 

The inverse matrix propagator in the leading order in the 
broken symmetry ground state is thus 



\{-iuJn + e-k) 

\{i(jJn+f-k) 

2\ Pc 




Here x = in B{K,T). The dispersion relation for uj is 
found by setting det = 0, which yields 



BiK,T) 



(29) 



After the analytical continuation iw„ — >■ w + iO+, the solution 
todetr>-i = Ois 



(30) 



Here A(w, k, T) = A/[l + \B{uj, k, T)] is the running cou- 
pling constant. We will show that B{uj,k,T ~ 0) — so 
X{uj, k,T = 0) = \. Therefore at T = 0, 



UJ' 



ek [ek + Xpc] ■ 



(31) 



One may compare this with the Bogoliubov dispersion o;^ = 
efc(efc + 2Apc) and see that the dispersions are similar. The 
factor of two comes from the fact that we have ignored the 
contribution from the anomalous density {(fxp) in the lowest 
order of our large-iV approximation. This contribution should 
get restored at higher order in the expansion. 

These results for the inverse propagator can also be dis- 
cussed in the 0(2) language by writing (j) = + i(f)2 and 



(/>* = 01 — Here we use the 0(2) symmetry to choose 
the condensate in the "1" direction. Then, in the broken sym- 
metry phase (j)^^^ — + a and 02^'' — n. The inverse 
propagator in the (x, cr, tt) representation takes now a slightly 
different form. The condensate density in this case only cou- 
ples to (J but not TT. We define ^ = (tr, tt, x)^ and the fluctu- 
ations are ^^V^^^. Using Eq. p?]) we obtain 




(32) 



Note that the time-derivative terms in Eq. p4| ) becomes 
TrdrCr — adrTT and this results in off-diagonal elements in the 
sub-matrix corresponding to tt and cr. From det 2?^^ = 
one finds exactly the same dispersion as the one given by 
Eq. ([30|. Therefore the Bogoliubov-like dispersion at T = 
emerges when one calculates the propagators in the correct 
broken symmetry ground state. 

Now we show B{K,T) explicitly. We define Gn = 
l/{—in^ + ujq) and G22 = + ooq), where ooq = 

eq + X. Then B{K) = {1/2)Y.q[Gii{Q)Gii{Q + K) + 
G22{Q)G22{Q + K)]. After summing over the Matsubara 
frequency, B{K, T) becomes 



IE 



l{ujq+k) - n{ujq) n{ujq) - n{ujq+k) 



(33) 



Here we have used n(a;g+ia;„) = n{u:q) and n(x)+n(— a;) — 
— 1. Defining Aw = u)q+k — and An = n{ujq+k) — n^uiq), 
we obtain the result 



b{k,t) = IJ2^' 



(34) 



The function B{uj, k, T) is then evaluated by the analytic con- 
tinuation iojn io + iO^ [33 1. One important consequence 
follows immediately. In the broken-symmetry phase LUq — eq 
and n{ujq) = as T ^- 0. Therefore B{uj, k,T = 0)^0 and 
this leads to a Bogliubov-like dispersion for the gapless mode 
inferred from the pole of the inverse propagator V^^. 

At finite T one can use the identity l/{x + iO"*") = 
P{l/x) — iTr6{x), where P denotes the Cauchy principle in- 
tegral, to obtain the full expression: 



B{u,k,T) = Jp^Ar 



(jj — Auj uj + Aw 



^ An[(5(a; + Aw) - 8(uj - Aw)]. (35) 



Thus at finite T one has to solve Eq. 
find the (complex) dispersion. 



with B(w, fc, T) to 
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III. THE NORMAL PHASE OF A TWO-COMPONENT 
BOSE GAS 

The effective action of a two-component Bose gas is a gen- 
eralization of Eqs. (|2]) and Q 



Al2((?!>t(?!>l)(02</'2) 



(36) 



We again introduce the large parameter N into the theory 
by the replication trick — ^ (/)j.„ where n — 1,2,...N 
and rescale the coupling constants \j — ^ Aj/A^ and A12 ^• 
X12/N. In the following we use a similar set of symbols to 
denote physical quantities of a two-component Bose gas. This 
set of symbols should not be confused with those for a single- 
component Bose gas in the previous discussion. 

The action with the source term after the replication be- 
comes 



S = / [dx 



N 



2^ "0 ^ 2N 



h.n + 



2N 



\n=l 




(37) 



Here we define $ 



5l,l,<?'l.l><?'2,l>'P2,ir 



= diag{h^+\h[\h^+\h[-\---). = 
Gq ^ — diag{^i, ^i, fi2, 1^-2, ■ ■ ■) is the bare (noninter- 
acting) Green's function of a two-component Bose gas, 

/i^*'' = zLdr — V^/ (2raj) foxj = 1, 2. There are N copies in 
(5(7^, and Gq ^. J is the source coupled to <i?. The identity 



1 = A/'M Vxi'Dx2DaiVa2 exp 



TV 
A^ 



— ! ^^2 



TV 



TV 



^2,n'P2,nJ 



(38) 



has the effect of introducing two delta functions similar to 
the case of a single-component Bose gas so one can replace 
Y^n=i^*j,n'i'3,n by {N/Xj)aj in S. This replacement facil- 
itates our resummation scheme and we will treat as a 



small parameter Let Gq 



Gn 



-diag{xi,Xi,X2,X2, ■ 



After integrating out 0j,n, one has 
1 TV 



S, 



eff 



dx] 



TV 
2X2 
1 



--J^GnJ- 
TVA12 



Ai 



TV 
A^ 



TV 



TV 



Ai2a2 + 
TV 



2AT"i^ 



*2 V^"l"2 - ^Xl"! 
A1A2 Ai 



-X2a2 



Tr In Gn^ - X 



(39) 



Here X = {xi, X2, cti, 0^2)"^ with its source term K and ex- 
pectation value Xc As in the single-component case we evalu- 
ate the path integrals over Xj , ctj via the method of stationary 
phase or steepest descent and in the leading order in large-TV 
we keep only the contributions at the stationary phase point. 

The generator of the one-particle irreducible diagrams is 
obtained from the Legendre transform of Seff- Explicitly, 
r = J{J^^c + K^Xc + Seff),vjheie^c is the classical value 
of <i>. We define the effective potential as Ve// = T/NVp. 
Keeping the leading term in the 1/TV expansion, and then set- 
ting TV = 1 we obtain the effective potential for static homo- 
geneous fields 



Ai' 



2Ai 



1 

2A^ 
1 



— — aia2 - T-XiQ!i - T-X2a2 
A1A2 Ai A2 



-Tr In Gn 



(40) 



Here $ 



)^ and Gn ^ has been reduced to a 



4x4 matrix. Again the Legendre transformation introduces 



the expectation values of 



and 



to r and Veff via 



J = Gq $ for the expectation values. The broken- symmetry 
condition is determined from the condition that we have found 
the true minimum of the effective potential: Weff/Scj)* — 0, 
which becomes Xj't'j = 0- In th^ normal phase (pj = while 
in the broken-symmetry phase Xj = 0. In the normal phase, 
the first term in Ve/ / is zero at the minimum of the potential 
(which occurs at (pj = 0). 

In the following we will focus on the normal phase of 
the mixture state and consider pi = P2 = Po ^nd mi = 
1712 — fn, where po is the density of a non-interacting single- 
component Bose gas with the BEC transition temperature 
To = 2Trfi^pl^^/[C/^{3/2)kBm]. Similar to the case of a 
single-component Bose gas, we define fco — and use fcg^^ 
and ksTo as the units of length and energy. 

The last term in Vef / can be evaluated using the standard 
Matsubara frequency summation technique and it becomes 
Y,k,j [wj72 + {1/(3) ln(l - e-^"^ )], where uj^ = ej + Xj and 
Cj = lT?k'^ / {2mj). To express V^jj as a functional of aj and 
Pj, we use 5Vef / /5aj = to obtain 



A12 

Xj = -^^J+aJ + ^aj- 



(41) 
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Here j — lif j = 2 and j = 2 if j = 1. This leads to 

3 ^ 12 



Al2 



-aiQ!2 



E 



ln(l 



(42) 



The renormalization of Vpff is similar to the proce- 
dure of a single-component Bose gas. Firstly one can 
show that f / SajSaj — —1/Xj+ (finite terms) and 

f / 5ai6a2 = — A12/A1A2+ (finite terms). This implies 
that the physical coupling constants only get finite renormal- 
ization, and as in the single field case A^ and A12 are equal 
to their renormalized values at T = 0. Then one may let 
\j — 4Trh?aj /rrij and A12 = 2iTh'^ai2/mr, where ai, a2, ai2 
are the s-wave scattering lengths of the intra- and inter-species 
collisions at T = and m.^ = niim2 / {mi + 1112) is the re- 
duced mass. To render the theory finite, we only need to con- 
sider the (infinite) renormalization of the chemical potential 
and vacuum energy. 

To make this procedure more transparent, we use Eq. ( |4T] i 
to express Vef / in terms of Xj for the moment. This gives 

Veff = Vo+J2 X3<l^3^3 + ^(Xi + Ml)' + 



A. 

2A2 

E 

k,3 



(X2 



- - ^(Xl + Ml)(X2 + M2) + 



1 



(43) 



Here = A^j — A1A2 andujj = ej+Xj- The renormalization 
of fij follows the set of equations 



A 



A2 A12 v;^ 1 A2 A12 

Ml - ^M2 + 2^2^ 1^^^^ ~ ^M2i? 

k 

Ai A12 1 _ Ai A12 

^M2-^Mi+2^2-^ 



^^2R - -TV^^^R■ (44) 



A2 



This renormalization absorbs the divergent term J^k jiXj/"^) 
in Veff - Then the vacuum energy is renormalized by 



Ai 2 A12 



2 ,,2 



Ai 



A12 



VoR + T^JJ-Ir + T^M2_R - TTMlflM2fl- 



A2' 



A2 



(45) 



Here we drop the subscript R and the vacuum energy. We 
first consider the normal phase of the mixture state. From 

6Veff/Saj = and pj = ~5Vef f / Sjij we obtain 



12 



a;"^ + a;a2"^ 



E 

k 



I \ I ^12 . 



(47) 



The solution along with 0j = then determines the extremum 
of Veff- To determine the stability of the mixture state, we 
compare the results with those obtained from the phase sep- 
arated state. Since our formalism uses the grand-canonical 
ensemble, one has to compare different states with the same 
chemical potential jij . The state with lower Vef / should be 
energetically stable. When the two curves of Veff intersect, it 
signals a phase transition into a different state. 

The broken-symmetry phase emerges when Xj vanishes ac- 
cording to the condition Xj4>j = 0. By analyzing ( |47] i with 
Xj = —lJ.j + oij + (Ai2/Aj)q;j one can see that this con- 
dition determines the critical temperature T™^ and for each 
component it coincides with the BEC transition temperature 
of a non-interacting Bose gas with the same density. For the 
case pi = p2 — Pa and nii — m2 — m, T™^ — T^™™ — Tq, 
which is independent of Ai , A2 , A12. One has to include higher 
order corrections in the large-A^ theory to get corrections to 
the transition temperature. 

The effective potential and equations of state for the phase- 
separated state are similar to those of the single-component 
Bose gas. For one of the species occupying part of the space, 
its effective potential is 



K7/ = (-M.+«r)(C)*c 



I/O, 

2A-("^- 



s\2 



(48) 



Here p,j needs to match the chemical potential of species j in 
the mixture phase. As a consequence, the density of the phase- 
separated state will be different from pj so we denote it by p^ ■ 
The energy dispersion is = ej — jij + a^- Since the BEC 

transition temperature scales as (p^^)^/'^, it is possible that 
in order to match pj, the phase-separated state may enter the 
broken-symmetry phase. Therefore we show the equations of 
state of the phase-separated state in the normal phase as well 
as in the broken-symmetry phase. 



In the normal phase. 



and 



This absorbs the divergent term J^k ji^j/'^) there is no di- 
vergence in Vef f after the renormalization. 

Following Eq. we let Xj — ^l^jR + ctjR + 

(Ai2/Aj)ajfl and rewrite Veff in terms of ajR. The renor- 
malized Vef f is 



V 



= E(- 



Pj 



A12 



A12 

aIX 



-Q!ia2 



k,3 



ln(l - 



(46) 



ps 

A, 



E"(-r)' 



ps 



k 



(49) 



Cj — Pj + ■ In the broken-symmetry phase. 



Here tjf - cj , 
p^j = {(jy^'')* (j)^" and one has 



ps 



ps ps 
P3 ^Pc,3 



k 



(50) 
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Figure 2. T4/ / as a function of A12 /Ai at T/Tq = 3 (evaluated at the 
mininum) for (a) feo^i = 0.04 and (b) fcoai = 0.1. Black (red) lines 
corresponding to the mixture state (phase-separated state denoted by 
PS). The dashed lines in the phase-separated state indicates that it is 
in the broken-symmetry phase. 



The dispersion is w^" = ej due to the broken-symmetry con- 
dition Xj — + o.j = 0. 

We now focus on the case where Ai = A2 = ATTh?ai/m. 
Figure |2] shows Veff from the mixture state and phase- 
separated state at T/Tq — 3 for two selected intra-species 
interaction strengths fcoOi = 0.04 and 0.1. For small A12/A1 
the mixture phase is more stable due to its lower Veff- As 
A12/A1 reaches a critical value, the two curves of Veff inter- 
sect and above the critical point the phase-separated state is 
more energetically stable. In the grand-canonical ensemble 
implemented here, the two states are compared at the same 
chemical potentials. Therefore the densities may not be the 
same in the two states. Note that when A12/A1 gets larger, 
the density in the phase-separated state increases in order to 
match the chemical potentials in the mixture state. Since there 
is no shift in Tc from the leading-order single-auxiliary-field 
theory when compared to a noninteracting Bose gas, the criti- 
cal temperature TP" = 2nh'^{pP'')^/^/[C^/^{3/2)kBm] of the 
phase-separated state increases accordingly. Eventually the 
phase-separated state may enter the broken-symmetry phase 
if T is not too high and we show this effect as the dashed lines 
inFig.|2] 

By locating the critical value of A12/A1 where the two 
curves of Vef / intersect at fixed T, we found the phase dia- 
gram shown in Figure [3] Each curve corresponds to the crit- 
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Figure 3. Phase diagram of the normal phase of a two-component 
Bose gas. The lines show the critical value of A12/A1 where a phase 
transition occurs for fcofli =0.1 (black) and 0.04 (red). Below the 
critical line the system is a mixture of normal gases and above the 
line a phase-separated state emerges. The dashed lines indicate that 
the phase- separated state is in a symmetry-broken phase. Inset: The 
normalized density difference {pps — p) / po at the critical line. 



ical line separating the mixture state and the phase-separated 
state. One can see that the mixture state prefers lower A12/A1 
and higher T while the phase-separated state prefers the op- 
posite. The dashed lines in Fig. [3] indicates that the phase- 
separated state enters the broken-symmetry phase. In that 
regime the leading order in large-iV approximation is only 
a qualitatively accurate result. In the region near and below 
Tc one can use the LOAF theory ifTTl ITSll to improve on the 
result presented here since that approximation exactly repro- 
duces Bogoliubov's results at weak coupling as well as cor- 
rectly includes the anomalous density and predicts a shift in 
Tc from the free-gas result. However, using LOAF will not 
change the answer in the region when both states are in the 
normal phase and would unnecessarily complicate the sim- 
plicity of the calculation presented here. One could also in- 
clude the next to leading order 1/A^ terms to be able to access 
the regime around Tc where the anomalous density correla- 
tions become important. A structural transition from a ho- 
mogeneous mixture state into a phase-separated state in the 
normal phase has also been studied in two-component Fermi 
gases with population imbalance [34 1- The underlying mech- 
anisms are different: For fermions the system is maximizing 
the pairing energy while for bosons the system is minimizing 
the repulsive interactions. 

Importantly, only the global translational symmetry is bro- 
ken in the phase transition from the mixture phase to the 
phase-separated phase when both species are in the normal 
phase. In the phase-separated phase, there is an interface sep- 
arating the two components and each component respects the 
local translational symmetry away from the interface. The dif- 
ferent densities of the two states across this mixture to phase- 
separation transition remind us of the liquid-vapor transi- 
tion, where no symmetry is broken and the density difference 
serves as the "order parameter" distinguishing the two phases. 
We thus study the normalized density difference (p^" ^ Pj)/ Po 
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at the critical line in the inset of Fig. |3] If the particle num- 
ber is conserved when one compares the two states, the nor- 
malized density difference should be 1 because the density in 
the phase-separated state should be twice as large as that of 
the mixture state. However, since we are comparing the two 
states at the same chemical potential, one can see from the 
inset of Fig. [3]that the conservation of the particle number is 
not respected. 

To draw the phase diagram with particle number conserva- 
tion, one has to work in the canonical ensemble with fixed 
particle numbers and find the corresponding free energy. The 
physics should be the same if the results are compared cor- 
rectly. For an isolated atomic cloud, our instability analysis 
may apply to a small region with the rest of the cloud treated 
as a reservoir. The phase separation could start growing if the 
mixture state is unstable in that focused region and the insta- 
bility may propagate to the whole cloud. 

One has seen that the instability of a mixture of two- 
component Bose gases can be analyzed using the mean-field 
approximation derived from the leading order in the large- 
N expansion, which involves the introduction of an auxil- 
iary field related to the normal density. The detailed struc- 
tures which develop when the system evolves into a phase- 
separated state, however, require numerical simulations of the 
equations resulting from the effective action and are beyond 



the scope of the present paper The width of the interface sep- 
arating the two species may be estimated using a variational 
method related to the one used in the estimation of the width 
of the interface separating two BEC phases in the ground state 
discussed in Ref. ifTTIl . 



IV. CONCLUSION 

We have shown that the leading order in our large- ap- 
proximation, which utilizes a single auxiliary field related to 
the normal density, leads to a mean-field theory usable at all 
couplings and temperatures that is a valuable tool for inves- 
tigating the physics of interacting Bose gases. For a single- 
component Bose gas we show that, by constructing the prop- 
agators in the broken symmetry vacuum, a Bogoliubov-like 
dispersion indeed emerges. For a two-component Bose gases, 
this approximation predicts a normal-phase structural phase 
transition between a mixture state and a phase-separated state. 
One possible application of two-component Bose gases is to 
simulate cosmological dynamics |35|. Our theory may help 
extend this application beyond the low-temperature regime. 

The authors acknowledge the support of the U. S. DOE 
through the LANL/LDRD Program. C. C. C. and F C. thank 
the hospitality of Santa Fe Institute. 
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